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Abstract In this work, we introduce a random field in view of natural image modeling, 
obtained as a limit of sequences of dead leaves models, when considering arbitrarily small or 
big objects. The dead leaves model, introduced by the Mathematical Morphology school, 
consists in the superposition of random closed sets, and enables to model the occlusion 
phenomena. When combined with specific sizes distributions for objects, they are known 
to provide adequate models for natural images. However this framework yields a small 
scales cutoff and a limit random field is introduced by letting this cutoff tend to zero. We 
first give a rigorous definition of the dead leaves model, and compute the probability that 
n compacts are included in distinct visible parts, which characterizes the model. Then, 
we derive our limit model and some of its property, and study its regularity. 

1 Introduction and motivations 

The structure of natural images is very specific, and strongly differs from the one of 
speech signals for instance. Most statistics of natural images exhibit non-gaussianity, as 
well as scaling properties. These two phenomena may for instance be easily observed on 
the distribution of the gradient of images gray levels ( |BB94| , |TMN98j ) . Other quantity 
bearing these properties include the power spectrum f |Kre52| . |Fie87| ) . wavelet coefficients 
( |Sim97j , |HM99j ) , morphological quantities ( |AGM99| ) or the distribution of local patches 
( |GK99j ). Non-gaussiannity is strongly related to the occlusion phenomenon. Indeed, in 
the process of image formation, objects hide themselves depending on where they lie with 
respect to the camera, which differs totally from an additive generation. This phenomenon 
leads to peculiar twodimensional structures such as homogeneous regions, borders and T- 
junctions. Besides, the scaling properties of an image may also be seen as a result of 
scaling properties present in nature. 

Several study f |Rud97| . |AGM99j . jEMHOU) show (either theoretically or experimen- 
tally) that most of natural images statistics may be reproduced though the use of a simple 
model of images, consisting in the sequential superposition of random objects, the dead 
leaves model of Mathematical Morphology. The mere nature of the model enables the re- 
production of characteristic structures of natural images (onedimensional discontinuities, 
homogeneous zones). Moreover, the use of a power law distribution x~ a for the sizes of 
objects enables to reproduce scaling phenomena (note that this distribution of object sizes 
was also considered in |Chi98j , but without occlusion) . In |LMH0l| , a version of the model 
corresponding to strict scale invariance is considered (that is a = 3) whereas |Rud97j and 
|AGM99j consider a as a parameter of the model. In both cases, the presence of small 
and large scale "cut-off" sizes is assumed, that is to say that objects sizes are bounded. 

In this article, we present a new model for natural images, that is obtained from a 
dead leaves model with scaling properties when letting the small scale " cut-off" frequency 
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tend to 0. By doing this, we model the small scales properties in a non-trivial way, and we 
are then in a position to study the regularity of images from a functional analysis point 
of view. 

In the first part, we recall some results on random closed sets and random tessellations, 
and then introduce the notion of a "colored" (or "textured") tessellation, and give a 
characterization of the convergence of a sequence of such random fields in the sense of finite 
dimensional distributions. Then, we define the dead leaves model, originally introduced by 
G. Matheron, [Mat68j . and study it in the framework of random tessellations. In particular 
we derive a new result enabling us to completely characterize the discontinuity set in the 
sense of random closed sets. In the last part, we introduce the specific dead leaves model 
with scaling properties, and derive some of its property in the case where objects are 
bounded. Then we look at the convergence of this model when objects sizes tends to 0. 
We first show that such a limit model is not well described by geometric means (i.e. by the 
theory of random closed sets) because there are small objects everywhere. Then, in some 
cases, we show that the finite dimensional distributions of the colored model still exhibit 
interesting behavior, and thus study its convergence from this point of view. Eventually, 
we study the regularity of this limit process using Besov spaces. 

2 Random closed sets and random tessellations 

In this section, we recall some facts about random closed sets, and give a definition of 
a random tessellation (that is a random covering of ]R 2 by disjoint sets) that enables to 
take into account covering by nonconnected sets. We then define the notion of a "colored" 
tessellation, and study the finitedimensional convergence of sequences of such processes. 

2.1 Random closed sets 

We first recall the definition of random closed sets (RACS). Let J 7 , Q and /C be respectively 
the sets of all closed, open and compact sets of W 1 . Let us set for any i C R", 

T A = {F G T : F n A = 0}, 
T A = {F G T : F n A / 0}. 

We define a topology Tjr on J 7 , generated by the basis of open sets 

In what follows, we write Bjr for the Borel <r-algebras generated by Tjr. 

Definition 1. A random closed set ofW 1 is a measurable function from a probability space 
(n,T,P) into (F,B T ). 
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We have the following characterization of RACS (see |Mat75j ^: 

Proposition 1. Bjr = a(J-Q,G G Q) = o~(J- K ,K G K). In particular, a random closed 
set F is thus characterized by its capacity function defined for every compact set K by 
G F (K) = Pr(Fr\K = 0). 

proof Using that for any K G /C, there exists a sequence (G n ) in Q such that G n J, K 
and that, for all G G Q, there exists a sequence (K n ) in K, such that K n j G, we obtain, 
respectively, for all K £ K, T K G <r(.F G , G G 0) and, for all G <E G, Tq G K G K). 

Hence the result. □ 

2.2 Tessellations 

A tessellation of 1R 2 is defined as follows. 

Definition 2. Let T = {Fi} be a collection of closed sets. We say that T is a tessellation 
ofR 2 if 

(i) U*i = R 2 . 

o o o 

(ii) for all i ^ j, Fi n Fj = 0, where F denotes the interior of F. 
(Hi) Card{i : Fi n A ^ 0} < +oo /or aZZ bounded icK 2 . 

Moreover we define the boundary ofT by dT = \^} i dFi, where dFi denotes the topological 
boundary of Fi . 

Observe that, by (jmj), the number of non-empty sets in a tessellation is necessarily 
countable. Hence, without loss of generality, we may index a tessellation by integers. Let 
T be the set of all tessellations indexed by N. 

Lemma 1. We have T G B^n, where B-pm is the a-field generated by the cylinders of the 
product space 

proof We notice that ((0) and (|ni))) is equivalent to 

n 

VN G N,3n G N, (j Fi D B(0,N), (1) 
i=l 

and 

ViV G N,3n G N,Vp G N, Card{i G {!,... ,p} : F f n J3(0, N) ^ 0} < n, 
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where 5(x, r) is the open ball centered at x with radius r. The second condition clearly 
define sets B-p®. Now we know from |Mat75| section 1-2] that (F, F') mFuF' is contin- 

o 

uous and thus is measurable and that F \— > (F) c is lower semi-continuous which implies 
measurability. Now note that condition (JnJ) is equivalent to 

FiHFj cdFiUdFj. 

We have ( |Mat75j ) that (F,F') — > F n F' is upper semicontinuous and that F — > dF is 
lower semicontinuous, and thus that these two applications are measurable. Eventually, it 
remains to prove that (F,F') — > l_F C pv is measurable. Note that F C F' is equivalent to 
FnF' = F, and that the diagonal set {(F, F),F E J 7 } is a closed set since T is Hausdorff. 
This proves the requested measurability. □ 

This enables us to define a random tessellation as follows: 

Definition 3. A random tessellation of the plane is a measurable map T from a proba- 
bilistic space (Q,S,P) to (3-^,13^), such that almost surely T G T. 

Property (iii) implies that for any random tessellation T, dT is a random closed set. 
In a random framework, there has been various definitions of random tessellations (see 
e.g. |Amb74j. |Cow80|. [MoTM]. |SKM9fij). mostly when the F { are convex polygons. A 
classical approach is to define dT directly as a random closed set without considering 
the underlying tessellation. However, it is not always possible to recover the i*Vs from 
dT (they may not be connected). Thus, we choose Definitional that provides a more 
informative definition of a random tessellation than its border. This is needed in Section 
12. 31 where we study the random function that arises when independently coloring each Fi, in 
a case where they are generally not connected. Classical examples of random tessellations 
(see the references in |SKM95| Chapter 10] and |OBSC00"| ) include Poisson line processes, 
Delaunay and Voronoi tessellations, and the dead leaves model that we consider in Section 

m 

2.3 A colored tessellation process and its finite-dimensional convergence 

In this section, we study the random process that arises when independently "coloring" 
(or texturing) each part of a random tessellation. In particular, we give simple criteria for 
the finite dimensional convergence of such a process. First observe that any tessellation 
T = {Fi} partitions M 2 into dT U (Ji(^A^ T )- 

Remark 1. Note that we could have given a more restrictive definition of a random tes- 
sellation where property (jmj) is replaced by 

for all i / j, Fi n Fj = 0. (2) 
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In this case, the plane is partitioned into dT U |J i i*V Also note that it is easily seen that 
if we assume that the sets F{ verify, 

for all i,Fi = F h (3) 

o o 

then (J2J is implied by property (indeed if x G Fj, there exists a ball B(x,r) C -Fj, so 

o o o 

that if Fi n Fj = 0, then x ^ -Fj). This gives an easy way to ensure that the boundary of 
the tessellation does not interfere with the interiors of the Fi, in which case some of the 
following definitions and results may be slightly simplified. 

Definition 4. Let T = {Fi} be a random tessellation and ({C(x) : x G K 2 }) be a 
real valued random field. For all x G M?\dT, denote by i(x) the unique index such that 

o 

x G F^y Let I = {/(x) : x G M. 2 } be the random field on IR 2 defined by 

f J(x) = C i(x) (x) /or a// x G M 2 \3T, 
\ J(x) = for all x G ST, 

where ({Cj(x) : x G M 2 }).; are independent copies of C, independent of T . We call I the 
colored tessellation process associated to T and C . 

The finite-dimensional distributions of / are easily obtained from the definition by 
conditioning on how the tessellation partitions the plane. The following notations allows 
for a formalization of this simple idea. For a given tessellation T = {Fi}, we denote by 

TZ^ the equivalence relationship defined by the partition IdT, Fi\dT^j , that is, for all 
x, y G M 2 xTZ^y if and only if either x,y G dT or there exists Fi such that x and y 

o 

both are in Fi\dT. For all n > 1, we denote by C n the set of all subsets of {1, . . . ,n} 
including the empty set. For any set A, we let 5(^4) denote the set of all partitions of A, 
with the convention 5(0) = 0. For all k G S(A) and all a G A, we denote by k(o) the 
equivalence class of a, that is k(o) = n(b) if and only if a and b belong to the same subset 
in the partition k. Let xi, . . . ,x n G M 2 . For all A G C n and for all k G S(A), we define 
the n-dimensional random variable 

I(A , K, Xx , . . . , X n J . — (C(A,K,l),...,C(A,K,n)), (4) 

where C(A, k, i) = if i <f_ A and C(A, n,i) = C K ^(xi) otherwise, {Cj, j G k(A)} being 
independent copies of C. Define the random variable (A, K) to be the element of the finite 
states space {(^4, k), A G C n , k G 5(^4)} defined by A := {i G {1, . . . , n} : xj ^ dT} and 
K is the partition of A defined by the equivalence relationship XjT^^Xj for all i,j G A. 
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Further define 
w(A, k, ( Xi )r =1 ) := Pr((A, K) = (A, «)) = 

(n n \ 

P| {i G A ^ x; £ 8T} n J { Xi ^ (T) Xj «(*) = . (5) 

i=l i,j=l / 

Lemma 2. For a// xi,...,x n G M. 2 , the distribution of (7(xi), . . . , I(x n )) is £/ie /i- 
ra£e mixture of {I (A, k, xi, . . . , x n ), ^4 G C re , k G 5(^4)} defined with respective weights 

{w(A,K,X 1} . . . ,X n ) 

,AGC n ,KG 5(A)}. 

proof Let xi, . . . , x n G M 2 . This lemma is easily obtained by observing that conditioning 
(J(xi), . . . , /(x n )) on (A, K) = (A, k) precisely give I (A, k, xi, . . . , x n ). Hence the result. 
□ 

In cases where any point has probability zero to belong to the set dT, the weights © 
are zero for A ^ {l,...,n}, which reduces the random vector (J(xi), . . . , J(x n )) to a 
mixture of random variables indexed by<S({l,...,n}). However, even in this simpler case, 
we will not pursue in studying the geometrical structure of the obtained distribution. A 
convergence result will simply be obtained by observing that the couple (A, K) defined 
in the proof above is a deterministic function (which we do not precise) of a random 
field indexed by (M 2 ) 2 which takes its value in {0, 1} (see Lemma 01 below). Let us first 
introduce the following definition. 

Definition 5. Let T = {Fi} be a random tessellation. For all n > 1 and for all 
xi, . . . , x n G M. 2 , let i?( n )(xi, . . . , x n ) denote the random variable which takes value one if 

o 

there exists i such that the n points xi, . . . ,x n are in Fi\dT and takes value zero other- 
wise. We will say that {i?( n )(xi, . . . ,x n ) : (xi, . . . ,x n ) G (M 2 ) n } is the rfi 1 order partition 
process. 

We obtain the following result. 

Lemma 3. Let Xi, ... ,x n G M. 2 and define the random variable (A,K) as in Lemma\^ 
Then 

(A^)ea(^'(xj) : x,y G{ Xl ,...,x n }). 

proof Clearly, (A, K) is a deterministic function of the partition processes taken at finite 
samples in xi, . . . ,x n , that is 

(A, K) G o-((R {m \z), z G {xi, . . . ,x n }"\ 1 < m < n}) . 
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Now observe that, for all m > 2 and for all yi, . . . , y n G M 2 , 

n-l 

fi (m) (yi,---,y n ) = n^ (2) (yi'yi+i)- 

t=l 

Furthermore, we have, for all y G ]R 2 , -R (1) (y) = # (2) (y>y)- The result follows. □ 

Thus it follows that the probabilities defined in © may be computed from the proba- 
bility distribution of a finite sample of the process . This allows for simple conditions 
to let a sequence (Ij)j£?q of colored tessellations converge to a limit field in the sense of 
finite-dimensional distributions. Let us recall that Ij — ^> Iqo if, for all n > 1 and for all 
xi, . . . x n G K 2 , (Jj(xi), . . . , 7j(x n )) converges to (/oo(xi), . . . , 7oo(x n )) in distribution. 

Theorem 1. Consider a collection of tessellations {Tj : j G N} and, for all x%, . . . ,x n G 
M. 2 , denote by R^ the second order process ofTj. For all j G N, let ({Cj-(x) : x G M 2 }) 
6e a sequence of real valued random fields, independent of the tessellations (Tj). Let us 
denote by Ij the colored tessellation process associated to Tj and Cj. Assume that 

(i) there exists a random process R® in {0, 1}( K2 ) 2 such that Bif ^ . 

(ii) there exists a random field Coo in M k2 such that Cj Coo . 

Then there exists a random field {/^(x) : x G M 2 } such that Ij — ^> loo- Furthermore the 

(2) 

finite- dimensional distributions of loo only depends on those of Rho and Coo. 

proof We use the expression of the finite distributions given in Lemma[2j Then assump- 
tion (0 and Lemma|21imply that, for all AdC n and for all k G S(A), Wj(A, k,xi, . . . ,x n ) 
(defined as in (JSJ) has a limit Woo{A, k, xi, . . . ,x n ) possibly computable from the finite 
distributions of i?oo . Assumption (JuJ) implies that, for all A G C n and for all n G S(A), 
Ij(A , K, Xx , . . . , X n ) (defined as in (gj)) converges to the random variable ioo(A «,xi, . . . , x n ) 
defined as in (J3J) using the random field Coo. Then (ij(xi), . . . , / 3 (x„)) converges in dis- 
tribution to the limit mixture {I^^A, k, x±, . . . , x n ), A G C re , k G S(A)} with weights 
{woo(A, k, xi, . . . , x n ), A G C n , k G S n }. Doing this for all n and for all xi, . . . , x n G M 2 , 

we see that all finite distributions converge and that the limit distribution only depend on 

(2) 

those of Roc and Coo- This gives the result. □ 

The two-dimensional distribution of the limit is detailed hereafter (as an application 
of the proof above). 

Proposition 2. Under the assumptions of Theorem^ the two-dimensional distributions 
of the limit process loo ar e given as follows. Let C'oo be an independent copy of Coo- 
For all x, y G M. 2 , (/^(x), ioo(y)) * s a mixture of the random variables (C 0C) (x), Coo(y)), 
(C 00 (x),C" 00 (y)), (Coo(x),0), (0,Coo(y)), (0,0) with respective weights Pr(i?^(x,y) = 
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1), Pr(B&\x,y) = O.i^fox) = l,R&\y,y) = 1), Pr(R®(x,x) = l,R®(y,y) = 0), 
Pr(i?^(x,x) = 0,i?^(y,y) = 1) and Pr(fi£ ) (x,x) = 0,R®(y,y) = 0). 

Remark 2. If, for all j G N, any point of M 2 has probability zero to belong to the set 
dTj, then, we have Pr(R^(x,x) = 0) = for all x G R 2 . Thus, for all x,y G M. 2 , 
(/(^(x), loo(y)) is a mixture of the two random variables (Coo(x), C^y)) and (Coo(x), C"oo(y)) 
with respective weights Pr(R^o (x, y) = 1), Pr(i?oo (x, y) = 0). 

3 The dead leaves model 

In this section, we introduce the dead leaves model. This model was introduced by G. 
Matheron in |Mat68j . and later presented by J. Serra in |Ser82j (see also |CT94j . j Jeu96j . 
|TK99Q . 

3.1 Definition 

Before proceeding with the definition of the model, we recall some notations that will be 
of use in the following. For any sets A and B, we set 

A = {-x,x G A}, 
AeB = {x£ M. 2 /x + B C A}, 
A® B = {x + y,x G A,y G B}, 

A 9 B is called the erosion of A by B, and A ® B the dilation of A by B. 

Since we will need these results in a near future, let us also mention that (F, K) — > 
F © K is a continuous (and thus measurable) function from T x /C into J- . We also 
have that (F, K) — > F Q K is upper semicontinuous from !F x K, \ {0} into J 7 . These 
results may be found in |Mat75j . Using the measurability of F — > from J- to and 
F ®K = (F c Q K) c , we see that (G,K) ^ G Q K and (G,K) ^ G® K are measurable 
maps from Q x IC into (Borelians are defined on Q in a way similar to those of F, 
see |Mat75| ). In the sequel, we will compute the probability of events involving such 
operations without any further comment. 

The dead leaves model is a particular instance of a random tessellation, obtained 
through sequential superposition of random objects falling on the plane. More formally, 
let {(xi,ti)} be a homogeneous Poisson point process (P.P.P.) in I 2 x 1 with intensity 
one. Let X be a random closed set of the plane, and (Xi) be i.i.d. closed sets with the 
same distribution as X, independent of the Poisson process above. Note that from the 
assumptions on {(xj,tj)} and (Xi), $ := {(xi,ti, X{)} is a P.P.P. of M 3 x F. We first define 
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Definition 6. The random closed set Xi + Xi is called a leaf and 




is called a visible part. 

Notice that, using the measurability of some standards operations on sets following 
from |Mat75j . {Vi} is a collection of random closed sets. From now on, we assume that X 
satisfies the following two conditions: 

(C-l) For all K £ JC, Ev(X © K) < +00, 

(C-2) there exists a disk D with strictly positive radius, such that Ev{X Q D) > 0. 

Proposition 3. M = {V{\ is a random tessellation o/IR 2 . 

In order to prove this result we will make use of the following two lemmas, the first of 
which will be repeatedly needed in the sequel. 

Lemma 4. Let K be a compact set, — 00 < t\ < t 2 < and define 

®K(ti,t 2 ) = {{xi,ti,Xi) : U C (h,t 2 ] and K C Xi + Xi} , 
<S> K {h,t 2 ) = {(xi,ti,Xi) : U C (ii,t 2 ] andKHxi + Xi ± 0}. 

Then $K(h,t 2 ) and ^ K (h,t 2 ) are P. P. P. 's of M 3 x T , respectively of associated mea- 
sures [iK,ti,t2 and /j, K ' tl ' t2 . Moreover we have /ix,ti,t 2 (^ 3 x = (*i — t 2 )Eu(XQK) and 
pKtuta^s x jr) = (t x - t 2 )Ev{X@K). 

proof The result follows by a classical result on independent thinning applied to <&, see 
|DV.T88j . and the fact that (K C (xi + Xi)) <^ (x { G Xj K), and (K n (x, + Xi) ^ 0) & 
{xi^Xi^K). □ 

Lemma 5. If K C M?, K measurable, is such that Ev[X Q K) > 0, then K is covered by 
some leaf Xi + Xi, for some i such that ti < 0, with probability 1. As a consequence, any 
bounded set is almost surely covered by a finite number of leaves. 

proof Let us fix t > 0. Using the same point process x as in Lemma HJ the probability 
that none of the leaves Xi with t < ti < satisfies K C + Xi is exp(— tEv{X K)), 
which yields the first assertion. Now let D be a disk such that Condition (C-EJ) is satisfied, 
that is Ev(X Q D) > 0. Since any bounded set K is covered by a finite number of disks 
with the same radius as D, the second assertion follows. □ 
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Proof of Proposition [3] We check that M is almost surely a tessellation. From Lemma 
we know that any compact set is almost surely covered by a finite number of leaves, 
which yields property (i) of Definition |2J Property (ii) of is directly obtained from 
Now, since any compact set K is almost surely covered by a finite number of leaves V% i , 
j = 1, . . . , n, letting T = min(i^ , j = 1, . . . , n), we see that all leaf Vj with tj < T has 
an empty intersection with K. It follows that the number of leaves intersecting K is less 
than a Poisson r.v. with intensity fi K,T '°(M. 3 x T) = T Ev{X@K) (see Lemma0} which is 
finite by Condition (C-HJ). □ 

o 

Remark 3. Recall that in Remark ^ we pointed out that if the Vj's are such that Vi = Vi, 

o 

then the Vi's verify relation (J2J. In the case of the dead leaves model, we easily see that 

o o 

if Xi = Xi, then the Vi's satisfy condition (j2J (even though in general Vi ^ Vi). We thus 
have a way to ensure that the interior of the Vi and the boundary of our tessellation are 

o 

disjoint, by ensuring that the X^s satisfy (J2J). Let us point out that the operation F — > F 
bears some similarity with the so-called "morphological opening", defined for e > by 
F — > (F Q B(0, e)) B(0, e), which is commonly used in image processing to "clean" sets 
or images. The aim is the same in our case : we can ensure that our objects do not have 
parts with empty interior, such as points or lines. 

We are now in a position to define a dead leaves model associated to a RACS X 
satisfying (C-Q) and (C-EJl as a random tessalation. 

Definition 7. (The distribution of) M = {Vi} is the dead leaves model associated with 
(the distribution of) X . The set dM is defined as in Definition^ and is called the boundary 
of M. The Vi 's are called the visible parts of M, and connected components ofM. 2 \dM are 
called the cells of M . 

Remark 4. We could also have defined the model through the random closed set dM. Here 
we propose a more complete definition since the tessellation M is generally not uniquely 
defined by its boundary (see the comments after Definitional. 

There exist different approaches for defining dead leaves models. The choice of intensity 
one for the P.P.P. {(xj,tj)} clearly is not restrictive. On the one hand, rescaling the x^s 
is equivalent, up to a global rescaling of the model, to a rescaling of X. On the other 
hand, any increasing transformation of the V s is unimportant as seen from the definition 
and, in particular, stopping the superposition of objects at time t = is not restrictive. 
Let us mention that M can also be defined as the stationary distribution of the Markov 
process (Mt)t>o, where Mt is the collection of visible parts associated to the thinned P.P.P. 
<f> R2 (-t, 0). Similarly, for a given compact K C M 2 , M f\K := {ViHK}, has the stationary 
distribution of the Markov chain (KnM t .^)i>Q, where (U(K))i>o is defined as the (almost 
surely) unique increasing sequence such that to(K) = and for all i > 0, ti + \{K) is the 
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Figure 1: Left : simulation of a dead leaves model, where the grain Xq is a disk with 
constant radius; Right : simulation of a dead leaves model, where the grain Xq is a disk 
with a uniformly distributed radius. 

closest tj strictly below U(K) such that (xj + Xj) n K ^ 0. Observe that K n M t is 
obtained by putting the new leaf Xj +Xj below the already fallen leaves. Hence the limit of 
(K n M t .t K \)i>o is hit at the first i such that K has been completely covered. It turns out 
that this i almost surely exists for all K (see Lemma EJ). The term "dead leaves model" 
originates from a more natural definition which consists in putting each new leaf above 
the previous ones and then consider the stationary distribution of this Markov chain. It 
is easily seen that the two definitions yield the same stationary distribution by a classical 
Markov Chain argument ( "coupling from the past") and that the former one, that is, 
the one hitting its limit, may be used for perfect simulation. This elegant argument was 
first introduced for the dead leaves model (and more generally for problems of stochastic 
geometry) in |TK99j , see also |Mol01| and the illustrating web applet |Kenj . The following 
simulations of the model have been computed this way, that is, by simulating the chain 

(K n M t . {K) )i> . 

In Figures ^ and|5J we show realizations of dead leaves models, for different choices of 
X. To visualize the model each grain is associated a color uniformly drawn in a table, and 
the model are "perfectly" simulated using the method described in the previous paragraph. 
In the first figure, X is a disk of constant radius (left) or a disk with uniformly distributed 
radius (right); in the second one, X is a rectangle whose sides are rocos(#) and rosin(0), 
with ro deterministic and 9 uniform between two bounds, uniformly rotated between 
and 7r (left) or a more complicated shape whose size follows a uniform distribution (right). 
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Figure 2: simulations of dead leaves models. Left: the grain Xq is a rectangle with a 
rotation uniformly distributed in [0,7r]. Right: the grain is more complicated shape, the 
distribution of its size is uniform. 



3.2 The functional Q 

The main practical result from the original paper by Matheron introducing the dead leaves 
model , |Matfi8| . concerns a functional, defined on the set of compact sets of the plane, 
equal to the probability that a given compact is included in a visible part of the model: 
for K G /C, 

Q(K) = Pr(3U <0,Kc V t ). 

Considering simple examples of possible K such as bipoints or segments leads to valuable 
geometric information on the model. 

Proposition 4. Let M be the dead leaves model associated with a random closed set X, 
then the probability of K € IC to be included in a visible part of M is 

' Ev(X<sK) 1 ' 

In what follows, we generalize this result by taking interest in the probability that 
n compact sets are included in distinct visible parts. In particular this gives us the 
probability that any compact consisting of a finite union of connected components hits 
the boundary of the dead leaves model. We define 

Q {n) {K x , ...,K n ) = Pr(3t h <---<t ln such that Ki C V h , . . . , K n C V in ). 
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Proposition 5. Let M be a dead leaves model associated with the random set X, and 
{Vi} its visible parts. Let us denote 

fP{k u ...,K n ) = Eu(x © k x ) Y[Eu((xe kj) n (x © K 3 -x) c ) , (8) 

and 

n 

GP(K u ...,K n ) = l[Ev(X®K j ), (9) 

i=i 

where, for all j, 

j 

Kj =\jK k . (10) 

k=l 

Then, the probability that n non-empty compact sets K\, . . . ,K n be sequentially included 
in n distinct visible parts is 



(11) 

G^(K 1 ,...,K n ) 

Remark 5. We actually see from the formula above that, because A Q B C A © C for 
all sets A, B, C such that B C C, we always have F^\Ki,. ..,K n )< G%\k 1} ...,K n ). 
Note also that (C-EJ) implies Eu(X) > and thus that G^p(Kx, . . . ,K n ) does not vanish 
for non-empty compact sets. 

proof We fix n compacts K\, . . . , K n . Summing over disjoint events we have that 

Q {n) {Ki, ...,K n ) = E\S2 < • • • < f[ M K j C V t 

where the sum is taken over all n-upplets of points in <!>. First note that from the definition 
of visible parts, the RHS of this equation may be writen as 

e i <■■■< tin) n c + ^)) n n + Xi 

\ j=l U>U. 

In the simplest case n = 1, this amounts to say that Q^ 1 ^ is the probability that there 

o 

exists a leave X, such that i^i is included in Xj and is not hit by subsequent leaves. We 
will now apply Campbell theorem to compute this expectation, and therefore need the 
following notations. Let £ := M? x R x T. We will denote by z = (x,t,X) the points of 
£. We write AA n ) (J\f for n = 1) for the space of point processes on £ n . For all n > 1, we 
then introduce the following point process on £ n , 

= {(^)" =1 : ^i, . . . ,z n e $ disjoints} . 
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We consider the following function, defined for all Z = (zj) T j =1 € £ n and for all $ = 

{(x lyj ,ti, v X iyj )] =1 }eM^ by 

i(h < ■ ■ ■ < i n ) f[ ( c + Xj)) n n (^J + = ) ) • ( i2 ) 

J=l \ *i,j6(tj,0] / 

In this equation and below, we use the convention that an indexed variable overlined with 
a tilda or a bar denotes a dummy variable as opposed to those related to When applied 
to ¥ = $W and = (^)"=i G $(n) » we obtain 

QM(K 1 ,...,K n )=El f(Z,^ n) ] 
Applying the refined Campbell theorem (see |DVJ88| ) . this expectation writes 



Q^(K 1 ,...,K n )= I f f(Z,^)P(dZ)V (z 1 ,...,2 n ,d^ 



where, writing Z = ((xj,tj, Xj))™ =1 , P(dZ) = YYJ =1 (dxjdtjP(dXj)), P being the proba- 
bility law of the random closed set X, and where Vq is the Palm distribution of the process 
<I>( n ) . Because is a P.P.P. in J\f, the Slivnyak's Theorem (see |DV.T88j ) applies, giving 

QW{K ll ...,K n )= ! I f(Z=(z j )] =1 ,($ + 6 Sl + --- + 6 Sn )W)P(dZ)P(d<!>), (13) 

J£" JAf 

where V is the law of process For all Z = (zj = (xj,tj, Xj))™ =1 such that t\ < ■ ■ ■ < t n , 
the product in (fT2|) may be written as 



f(z,(<s>+s Sl +---+5 Sn )W)= mi^c^ + i^l \^pi(K j _ 1 n(x j + x j ) = \ 

II II HK j n(x i + x i ) = ®)\ 11 i(K n n(x k + x k ) = $), (u) 

K j=itie(tj,t j+1 ] J t k e(t n ,o] 

with K_j as defined in (|10[). To compute the integral in (|13j) . we may first integrate the 
term of the second line of 1)14(1 . first with respect to V, then to l(ti < ■ • • < t n )dt\, . . . , dt n . 
Now, at fixed s < t < 0, and for K compact, we have that 

V{K r){xi + Xi) = ® for all U G (s, t]) = exp ((s - t)Ev{X © if)) , 
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which follows from Lemma ^ Then using a change of variable uj = tj — ij+i, for j 
l...n — 1, we obtain 



Q( n \K 1 ,...,K n ) = f[Ev(X®K j ) 1 



/ 



u=i / \i=2 / 3=1 



(15) 



The first term of the RHS of the previous equation is (G^ ) 1 , and the term of the second 
line writes 

TT ( f l(Kj C x + l)!^,! n (x + X) = 0)dxP(dX)^) , 
with the convention K_q = 0. Now, for two compact sets A and B, we have 

J l(4c(i + X))i(fi n (.x + x) = Q)dxP(dX) = e x v{{x ei)n(i® Bf), 

which, along with the last equations, yields pijl. □ 

In case n = 1, we get the original result of Matheron, Proposition^! and the case n = 2 
was treated in |.Teu96j . Note also that we can compute the probability for n connected 
compact sets K\, . . . , K n to avoid the boundary of the dead leaves model M. In case n = 2 
for instance, this is 

Pr((#i U K 2 ) n dM = 0) = QW(K U K 2 ) + Q (2) (K 2 ,i^) + g (1) (^i U K 2 ). 

More generally, we show that the knowledge of for all n uniquely determines the 
probability for any compact set K to avoid the boundary of M, which characterizes its 
distribution in (J 7 , Bp). 

Proposition 6. Let M be a dead leaves model, then the distribution of the random closed 
set dM is uniquely determined by the junctionals Q^ n \ n € N. 

proof Thanks to Proposition^ it is enough to prove that the functional Q(K) = Pr(KCl 
dM = 0), defined for compact sets K, is uniquely determined by the Q^ n \ Let K G /C, let 
r n > be a sequence converging to 0, and for each n, let {x\ n ^}i=i t ,,, t N„ be finite sequences 
in K such that K C C n = UiD(xf,r n ), where D(x,r) is the (closed) disk centered at x 
with radius r. Note that since each C n is a finite union of connected compact sets, the 
knowledge of the i £ N, uniquely determines Q{C n ). Now since C n [ K, we have 
that F Cn j T K , and thus that Q{C n ) | Q(K). □ 
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4 A dead leaves model with scaling properties 



In this chapter, we take interest in a dead leaves model that has a specific object distribu- 
tion. X, the grain of the model, is equal to RY, the homothetic of a random compact set 
Y, where R is a random variable with density /, independent of Y. We will consider in 
details the case f(r) = Cr~ a \{r$ < r < r±), with a > 1, that is when the size distribution 
of the objects satisfies some scaling properties within in a given range. As explained in 
the introduction, this choice for / is motivated by natural images modeling. Clearly we 
cannot take vq = for / to be a density. Now, from a modeling point of view, taking 
ro > is not satisfying as well because this would be inadequate for capturing the pres- 
ence of small objects in natural images. From a theoretical point of view, this reduces the 
model to only very simple smoothness classes (namely, piece-wise constant images). From 
a practical point of view, it means that there exists a minimal size for the objects in the 
image. It is not clear at all what physical meaning to give to this minimum object size, 
and how to deal with this supplementary parameter of the model. It is also unclear how 
to relate this minimum size to the resolution of an image, which we may assume to be 
obtained from the model through filtering and subsampling. Moreover, this contradicts 
empirical experiments (see IGMOlp which conclude to the presence of small "objects" up 
to the smallest observable scales in digital images. Therefore it is worthwhile to wonder 
about the limit of the model as ro tends to zero. The parameter r\ is not crucial for 
modeling smoothness properties because it does not influence the small scales behavior, 
except perhaps when r\ = oo or n tends to oo in which case the model may degenerate. 
These cases will also be discussed. In the good cases we will consider both the convergence 
of the model boundary set and the finite-dimensional convergence of a colored version of 
the model, giving raise to quite different kinds of limits. 

4.1 Definitions 

Let 7 C t 2 be a random compact set, and, for any < ro < r±, define the probability 
density function 



with r](ro,ri) = (1 — a) _1 (rQ~ a — r\~ a ), / = otherwise. We consider M(ro,ri) the 
dead leaves model associated with the random closed set X = RY, where R is a random 
variable with density /(ro, r±, •). For convenience, our notations do not refer to the scaling 
parameter a. However, it must be kept in mind that these definitions highly depend on 
this parameter. In order for this model to make sense, we want to keep Y within fixed 
proportions. All along the paper it is assumed that 

(A-l) there exist a-i > «i > such that almost surely D(ai) C Y C Dfa), 




(16) 
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where D(a) is the disk of radius a centered at the origin. 

It is straightforward to check that X defined this way satisfies the requirements for 
M(ro,ri) to be a dead leaves model, (C-0 and (C-|2J) for all < ro < n < oo and, if 
a > 3, also for n = oo. 

Let p(ro,ri,x) denote the probability that the origin and x 6 M 2 are in the same 
visible part of the dead leaves model M(rQ,r\). In the next section we focus on the 
limit behavior of p(ro,ri,-) as ro and r\ goes to and oo respectively. This will tell us 
under which condition for a a non-degenerate limit process can be obtained. According 
to Formula 0, for all x £ M 2 , 

Eu [XQ{0,x} ) Eu[Xr)(x + X 



pCro.ri.x)- £i/(Jf0{Ox}) - Ev ( X \j(x + X)) ' (17) 

where v is the 2-dimensional Lebesgue measure, and E the expectation with respect to the 
law of X. Writing Ey for the expectation with respect to the law of Y, Fubini's Theorem 
and then a geometric argument along with (|16|) give 



&(lnx + l) = | E Y v i^uY n (x + uY)j /(r , n,u) du 



v(ro,n) / 7(- J u 2 ~ a du, (18) 



r 

where 7 denotes the geometric covariogram of Y, see |Mat75j . that is, for all y £ 111 



7 (y) := ( y n (y + y) 

Let us now consider some basic assumptions and results concerning the mapping 7. From 
(A-^), we obtain 

Lemma 6. The mapping y 1— ► 7(y) is continuous over M 2 , 7(0) > 7ra 2 and, for ally € M 2 , 
< 7(y) < 7(0)1 where \ ■ \ is the Euclidean norm. Moreover, 

|y| > 2a 2 => 7 (y) = 0. (19) 

(o o \ 

proof The bounds on 7 are immediate. We now notice that y 1— > v I Y n (y + Y) I is the 



convolution of the indicator function on y with itself. Since y is bounded, this indicator 
function is square integrable with respect to v and the convolution is continuous. By the 
dominated convergence theorem, the continuity is preserved after taking the expectation. 
□ 

The two following lemmas show that a stronger regularity of the geometric covariogram 
at the origin may be obtained by imposing some geometrical properties on Y . First it is 
known, see |Mat751 Proposition 4-3-1], that 
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Lemma 7. Let u be a unitary vector. Let C be a closed convex set such that C ^ 0, then 
x i— ► v(C fl C xu ) admits a derivative to the right at x = equal to vi(P u (C)), where v\ is 
the one dimensional Lebesgue measure, and P u the orthogonal projection in the direction 
ofu. 

This lemma can for instance be generalized to finite unions of convex sets. Simple 
criteria on the random set Y can thus be deduced in ordered to impose that 

(A-2) for every unitary vector u, the function 7(0:11) is differentiate to the right at x = 
and the derivative is bounded away from the origin and infinity independently of u. 

This condition will be used as an assumption in some of the results below. However, 
in most cases we will relax this assumption into the following one. 

(A-3) For any 6 > 0, we have 7(x) = 7(0) + o(|x| 1_<5 ) when x — ► 0. 

For obtaining such conditions, we have the following bound. 
Lemma 8. Let C be a compact set. Then for all x € R 2 , 

< v{C) - v{C n (x + C)) < v I [J (u+[0,x])j , 

\uedc ) 

where [0, x] denotes the segment of points between x and the origin, that is {ax, a € [0, 1]}. 

proof Let x € M 2 . We may write 

< u(C) - v(C n (x + C)) = v(C\(x + C)) 

00 o o 

Let y € C\(x + C), that is y € C and y — x ^ C. Thus the line segment [y,y — x] 
intersects dC. The result follows. □ 

o 

In particular we see that if dC is a rectifiable curve with length L, then < v{C) — 

o o 

u{C n (x + C)) < L|x|. More generally, if C is such that its topological boundary, dC has 
an upper box-counting dimension less than or equal to one, then (sec Fal90 , Proposition 
3.2]), for all 5 > 0, 




where D(x,a) denotes the disk of radius a centered at x. By Lemma El the same bound 

OO o 

applies to v(C) — v(C fl (x + C)) with a = |x| in a case where the boundary of C is not 
necessarily a parameterized rectifiable curve. Hence Assumption (A-EJ) is easy to verify 
for a large class of random closed sets Y if one has some minimal property on the measure 
of dY . This is summarized in the following corollary. 
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Corollary 1. Assume that for any 5 > 0, we have 

E Y y(\J (u+[0,x])j =o(|x| 1 - 5 ) 
\uedY / 

when x — > 0. Then (A\$j) holds true. 

However, in order to illustrate that some smoothness is needed on the boundary of 
Y to insure (A-J5J), even for a deterministic Y. Consider the following example. Let 
h : [0, 1] — » (0, 1) be a continuous mapping and let us define 

C ■= {(x, y):0<x<l, h{x) - l<y< h(x)} . 

Then we have, for all u > 0, 

v(C)-v(Cn({u,0) + C)) = v(C\{(u,0)+C)) = u+ f \h(x)-h(x-u)\dx. 

J u 

Hence in this case, taking say Y = C non-random (it could be made random by taking 
h random), (A-OJ) would imply that, for all positive 5, \h(x) — h(x — u)\ dx = o(u 1_l5 ) 
as u I 0. This exactly means that h £ n«5<i 5 i ([0, 1]) (see e.g. |DL91] or |Mey90| and 



Section 1131 where the Besov spaces Bp' q are defined) and continuous bounded mappings 
may be easily found out of ^'^([0, 1]) for some 5 < 1. In such cases (A-I2J) cannot 
be satisfied. This example shows that this assumption amounts to some smoothness 
assumption on the boundary of Y. 

Another way to interpret Assumption (A-EJ) is to say that in the limit model we shall 
consider, the smoothness is mainly driven by objects sizes distribution and not by the 
irregularity that may occur in their boundaries. 

4.2 Basic convergence results 

In order to simplify the presentation of this section, we temporarily assume that the 
distribution of Y is isotropic. However, all the convergence results of this section remains 
true in the case of nonisotropic Y, the various quantities under study then depending on 
a directional parameter. In the isotropic case, we let 7 and p(ro,n, •) be mappings of the 
real variable x = |x|. The general case may be obtained in the same way by letting these 
mapping also depend on the angle coordinate of x. 

From (|18j) we may investigate under which conditions p(ro, r%, ■) has a non-degenerate 
limit when one wants to push the model towards the values of ro and t\ which are not 
allowed, that is, for all a > 1, ro = because /(0, r±,-) is not well defined and, for all 
1 < a < 3, r\ = 00, in which case, condition (C-^) does not hold. 

The cases for which the obtained limits degenerate are summarized in the following 
Proposition. 
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Proposition 7. For all x > 0, 



(i) if 1 < a < 3, lim inf p(ro, r\, x) = 1, 

ri— »oo 0<r <ri 

(ii) if a > 3, lim sup p(ro, n, x) = 0, 

ro->0 o<r <n 



(iii) if a = 3, lim 

r ->0 
ri^oo 



p(r ,ri,x) 1-2 



'log(ri) 



proof Using (JTBJ) and since ^(IU(x + I)) = 2i/(X)-i/(lnx + I), Eq. (T7J) reads 



p(r ,ri,x) 



/,?[7(0) + (7(x/u)-7(0))]« 2 - a d« 



(20) 



We now derive the asymptotic behavior of the integrals in this formula depending on the 
value of a > 1 . 

Take 1 < a < 3. Let x > 0. From Q19JI and since 7(0) > 0, for any e > 0, there exists 
uq such that for all u G [no, 00), 

| 7 (x/n) - 7(0)| <e 7 (0). 

Since the integral below diverges as r\ — > 00 and is bounded as ro — > 0, for n sufficiently 
large and for all ro < r±, 



[ ri u 2 - a du< (1 + e) / 



'1 
«oVr 



n 2 ~ Q du. 



From the two last equations and from (|20[). we get, for all ri sufficiently large and for all 

< r < n , 



1 > p(ro,n,x) > 



ma - e) Q vr y- a du 1 



7(0)(l + e) 2 J£ Vn> ua-«du (1 + e) 2 ' 



Hence (0) by letting e decrease to zero. 

Take now a > 3. From (fl9|) and (|20|) . we have, for all ro < r\ < x/(2a-z), p(rQ,r\,x) = 
and, for all ro < x/{2a,2) < r±, 



n 



j(x/u)u a du 



ro 



n 



x/(2a 2 ) 



[j(x/u)]u 2 - a du. 



Since 7 is bounded by Txa 2 from above, we get, for all ro < x/(2a 2 ) and for all r\ > ro, 



p(r ,r!,x) < 



27(0) /, 



x /i 2a 2) „,2 



ro 



II- 



du-va 2 f^ (2a2) u 2 a du 
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which does not depend on r\ and tend to zero as ro — > 0. This gives (JnJ). 

We now conclude with the case a = 3. From (|19|) and the continuity of 7, the numerator 
of the RHS of (|2U|) behaves as 7(0)log(ri) when ro and r\ respectively tend to and 00. 
For the same reasons, the denominator behaves as 7(0)(log(ri) — 21og(ro)). We finally 
obtain (jmj. □ 

It is worth commenting on the consequence of these simple convergence results. In 
case (0), the result says that, as r\ tends to infinity, however ro < r\ may behave, any two 
points end up in the same visible part. The big objects predominate at the limit so that 
the image is simply covered by one single set X with probability one. In case ((HI), the 
result is the exact opposite. As ro tend to zero, however r\ € (ro,oo] may behave, any 
two points never belong to the same object. The small objects predominate so that every 
visible part reduces to a point and the limit distribution of the image is white noise. See 
Figure 01 for an illustration of these cases. Finally, in the case (JmJ), the limit depend on 
the behavior of log(ro)/ log(ri). Convergence to 1 or as in cases (Q) and (jnj) are observed 
if only one of the limit ro — > or r\ — > 00 is taken. If, for instance we take ro = r[~ s for a 
fixed s, and let r\ tend to 00, we obtain a limit which depend on s but does not depend 
on x. 

We will generally avoid these cases in the sequel as they obviously give degenerate 
limits. In contrast, when ro — ► for 1 < a < 3, the limit of p is easily obtained from (|2()j) 
giving raise to the continuous prolongation 

Jo [ 2 7(0) - lix/u)] u 2 du 

These integrals trivially converge since 7 is bounded and a < 3. In the following result, by 
assuming sufficient smoothness on the boundary of Y, we exhibit asymptotical behaviors 
in which the geometry of the model only appear in multiplicative constants while the 
qualitative behavior is a power law of the form x 3 ~ a . 

Proposition 8. We have the following asymptotic expansions. 

(i) For all a > 3, p(rQ,oo,x) = g(a) (x/ro) 3 ~ a (l + o(l)) as x/tq — ► oo ; where 



a — 3 



CO 



g(a) := ^-(q) I ^(V^) a dv < 00 for all a > 3. 





(ii) Under (A\$, for all 2 < a < 3, l-p(0,n,x) = g(a) {x/rif~ a {l+o{l)) asx/n -> 0, 
where 

g( a ) ■= 2 ^~ a ^ f (^(0) _ 7(l/u)) v 2 ~ a dv < 00 for all 2 < a < 3 
7(0) Jo 

(iii) Under (A\E), for a = 2, 1— p(0, r\,x) = %^r- (x/r\) log(x/ri)(l+o(l)) as xlr\ — > ; 

7(0) 
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(iv) Under (A\^j, for all 1 < a < 2, 1 — p(0, r\,x) = g(a){x/r\)(l + o(l)) as x/r\ — > ; 
where 

g(a) := ^- — < qq f or a // i < a < 2. 
(a -2) 7 (0) 

In the formula above, under (A\$) and having assumed that Y has an isotropic distribution, 
7(0) denotes the right-sided derivative of x 1— > 7(x), x £ [0, 00), at the origin. 

proof Take a > 3. An obvious change of variable gives 

roc 

^{x/u)u 2 ~ a du = x 3 ~ a / j(l/v)v 2 - a dv. 

ro J ro/x 

Using (fl9|) and (|2()|). letting x decrease to zero gives (0). We now take a < 3. We similarly 
have 

ri pn/x 

( 7 (x/u) - 7(0))u 2 " a (in = x 3 " Q / (t(1/u) - 7(0))v 2 ~ Q dv. 
■/ 

From (|2*U)) and (|2*T|) and standard computations, we obtain ([n |) -((iv jl . □ 

It is worth noting that empirical experiments have shown that natural images exhibit 
such power law behaviors, with a having values usually in the range [2.5, 3] (see |Rud97j . 
and also various study on the power spectrum of natural images) . Hence in the following we 
will focus on cases (0) and (JuJ). One way to interpret these two cases is the following. The 
behavior at large scale of p(ro, 00, x) with ro > and a > 3, is qualitatively similar (namely 
a power law) to the behavior at small scales of p(0, r%, x) with r\ > and 2 < a < 3. Hence 
the limit of M(rQ,r\) as ro — ► for a fixed r\ > seems the natural way to extrapolate 
the model M(ro, r±) from ro > 0, r\ = 00 and a > 3 to ro = 0, r\ > and a < 3. Based 
on this observation, we dedicate the next sections to defining a limit model M(0,ri) for 
2 < a < 3. It turns out this limit is degenerate from a random set point of view (see 
Section [2.1(1 but it is not from a colored model point of view (see Section [2.3(1 . 

4.3 Limit of the boundary set 

In this section we assume that 1 < a < 3 and n > or that a > 3 and n £ (0, 00]. We 
have seen in Section 14.21 that, although the model M(0, r\) is not correctly defined, the 
probability p(ro,n,x) can be continuously prolonged at ro = for all x. In this section, 
we investigate how the model M(ro,n) converges under the same conditions. Here we 
take the classical point of view of random closed sets and we consider the distribution of 
the (random) boundary set. In the case a > 3, p(0, ri,x) = 1 (see Proposition EJ) so that 
we expect a degenerate limit of the dead leaves model. It turns out that the limit is also 
degenerate in the case 1 < a < 3 as shown by the following proposition. 
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Proposition 9. Let dM{rQ,r\) denote the random set consisting of the union of the 
boundaries of visible parts in the model M{rQ,r\). Then 

dM -» R 2 , 
r -»0 

the convergence being the weak convergence of random closed sets. 

This convergence result follows easily from the next proposition which investigates the 
presence of constant areas in a possible limit process as ro — > at fixed r\. 

Proposition 10. Let Q{r$, ri,r) denote the probability for a disk of radius r to be included 
in a visible part of M(ro,n). Then, for any r > 0, lim Q(ro, n, r) = 0. 

r ^0 

proof According to formula and then to (A-^), we have 

_ EviX Q D(r)) Ev(RD(a 2 ) Q D{r)) _ Jgv < ua * ~ r ) 2 u ~" du 
Q( - r ° ,n,r) ~ Ev(X ® D(r)) ~ Eu(RD{ ai ) © D(r)) ~ Q 7r(u ai + r) 2 u~ a du ' 

The limit is now obvious. □ 



proof of Proposition |9] Let P{rQ,r\, •) denote the probability law of dM(ro,r\) in the 
probability space (J 7 , Bjr) (see Section I2TT]) . We recall that a sequence P n weakly converges 

to P in (F,B?) if for all E G Bp such that P{E) = P(E), P n (E) converges to P(E) (see 
|Bil68j ). Moreover, in the case of the probability space (J 7 , B^), this amounts to check 
that for all K G K, such that P^k) = P(F<> ), Pni^x) converges to P{!Fk) (see |Lya83 , 

K 

|Mol93j ) . Here the limit distribution P associated with the deterministic set M 2 satisfies 
P{3~k) — 1 f° r an compact set K ^ and P(T§) = 0. Take a compact set K such that 

o 

K ^ 0. There exist a disk with positive radius r included in K so that 

P(ro,n, Tk) > P{ro,ri,FD(r)) = 1 ~ Q(ro,n,r). 
The result then follows from Proposition 1101 □ 

4.4 The colored dead leaves process and its limit 

For modeling the effect of the occlusions on the smoothness in natural images in a non- 
trivial way, we consider the limit of the model M(ro,ri) as ro — ► 0. From the point of 
view of random closed sets, the limit degenerates as we have seen in Proposition EJJ The 
possible limit set is not well described by the boundary limit, and we now take interest in 
the limit of a colored dead leave model in the more traditional sense of finite-dimensional 
distributions. Since a dead leaves model M is a random tessellation (composed by its 
visible parts (Vj), see Proposition |5J), we may use the results of Section [231 
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Definition 8. Let {Cb(x), x G IR 2 } be a random field. We denote by /(Co,ro,ri) the 
colored dead leaves model obtained from the random tessellation M(ro,ri) (see Definition 
0). If Co is the constant random filed with uniform marginals, that is, for all Xi, . . . ,x n 
and for all c±, . . . , c n € (0, l) n 

P(C (xi) < ci, . . . ,C (x n ) < Cn) = P(C (xi) < min(ci, . . . ,c n )) = min(ci, . . . ,c n ), 

we simply denote the colored dead leaves model by I(ro,r\). In other words I(rQ,ri) is 
obtained from the dead leaves model by independently coloring each leaf with a uniform 
distribution. 

Remark 6. Note that the definition of I requires the knowledge of (Vi), and not only the 
distribution of dM, since a visible part is not necessarily connected. 

Remark 7. Observe that, if Co is a stationary process, then the same is true for /(Co, ro, r±). 
For instance, I(ro,ri) is stationary. 

We now investigate the existence of a continuous prolongation of /(Co, ro,r{) at ro = 0. 
As in Section l4.3( we assume that 1 < a < 3 and n > or that a > 3 and n € (0, oo]. 
Recall that in that case the model M(0,r\) is not correctly defined although M(ro, ri) is 
for all ro S (0,ri). We have already taken interest in the limit of this model as ro — > 
in the sense of limit distribution of its boundary set (see Section 14.3ft giving raise in 
degenerate limits even when the first order probability function continuous prolongation 
p(0, r±, •) is non-trivial. We now investigate the limit a colored dead leaves process defined 
from M(ro,ri) and a given color random field Co- From Theorem^ for all random field 
Co, it is sufficient to study the limit of the corresponding second order partition process 
which we denote by {i?( 2 )(ro,ri,x, y) : xjER 2 }. The following remark and assumption 
will highly simplify this problem. 

Remark 8. Applying Proposition 01 we have, for every x S M 2 , Pr(x £ dM) = 1 — 

o 

Q({x}) = 1 — Eu((X))/E(iy(X)). This is clearly zero if we assume 

(A-4) u(8Y) = almost surely. 

This is always true for non- degenerate cases since dY typically is a curve. Note for instance 
that (A-|3J holds under the assumption of Corollary^ In the sequel we assume (A-0J) for 
convenience. 

We obtain the following result. 
Proposition 11. There exists a random process /(Co,0,ri) such that 

I(C , r , n) — > /(C , 0, n) as r -> 0. 
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proof Let n > 0. In this proof, for all ro G (0, ri), we denote by {V^ r °}i the visible 
parts of the dead leaves model M(ro,ri). Using Theorem ^ it is enough to prove that 
there exists a random process R( 2 \0,n) such that i?( 2 )(r ,n) ^i?( 2 )(0 , n) as r — > 0. 
More generally we may show that, for all n £ N, for all compact sets K\, . . . , -ff n , the joint 
distribution of ( l(3ij, C Vj°) ) G {0, 1}™ converges as ro tend to the origin. Then 



' 3=1 

applying this to compact sets composed of two points gives the result. Let us denote in 
this proof section, for all n > 1 and for all ei, . . . , e n £ {0, l} n , 

pW( £l> . . . , e n ) := Pr((l(3i i9 If, C l£>)) J =j = (e,)J =1 . 

Observe now that, for all n > 2 and for all ei, . . . , e n _i G {0, 

pW( eij . . . , en-1,0) + p (n) (ei, • • • , £n-i, 1) = P (n) (ei, . . . , e n _i). 

A simple induction on n thus gives that it is sufficient to show that the probability 
p( n \l,...,l) = Pr(3h,...,%n/Ki C V[°,...,K n C V[°) converges in [0,1] as r tends 
to 0. We may also assume without loss of generality that each Kj contains at least two 
distinct points. Otherwise Kj is included in the interior of a visible part with probability 
one (see Remark EJ). Finally, we claim that it is now enough to prove the convergence of 

o o 

Pr(3h, i n /K x C V[°, ...,K n C V[° and t h < t h < ■ ■ ■ < t in ). 

This follows from the fact that the union of two compact sets is a compact set, so that we 
may restrict ourselves to disjoint visible parts, by an elementary induction. The probability 
above is defined as (K±, . . . ,K n ) as in After Proposition [SJ we know that 

Q%>{K U ...,K n ) = FW(K U K n )/G%>{K U ...,K n ), 

with obvious notations. Pick a Kj and let 5 denote its diameter. Recall that we have 
assumed 5 > 0. From (A-Q), we have v{RY Q Kj) = for all R such that 2Ra2 is larger 
than 5. This implies that, 



v(ro,n 



)Eu {{X Q Kj) D{X(B Kj^f^j = jf 1 E Y v {{rY Q Kj) n (rY © Kj^f^j r~ a dr 



stays constant as soon as tq goes below <5/(2a2). Here, for j > 1, K_j is defined in (|10|) . 
the case j = 1 being obtained with the convention K_ = 0. Hence, from (JSJ), it is clear 
that, for ro small enough, {K±, . . . , K n )rj(ro, r\)~ n does not depend on ro, with rj as 
in (|16j) . On the other hand, for all j = 1, . . . , n, we have 

Ev{X © Kj) = rj(r ,n) / r~ a E Y v (rY © Kj) dr, 

Jr 
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where we recall that Ey is the expectation with respect to Y . Since the integrand is 
positive, 77(7-0, ri) _n Gro^ {K\, . . . , K n ) has a limit in (0, 00] (it is non zero since the /Q's 
are non empty). Simplifying by 77(7-0, n) - " in the ratio defining Q^ 1 above, we obtain 
that it has a limit as tq tend to the origin, which, as we claimed, is sufficient for showing 
Proposition [HI □ 

Note that in this proposition, we did not make any assumption on the value of a. 
However, we will see in Proposition El that this result is interesting for 1 < a < 3, in 
which case there exist a measurable version of our limit process. If a > 3, the limit process 
simply is white noise. 

We conclude this section by a simple corollary of Proposition 111! where we compute 
the two-dimensional distributions of /(Co, 0, n). We have, for all x, y e M 2 , 

Pr(i? (2) (r ,ri,x,y) = 1) = p(r , n, y - x), 

where p(tq, n, •) is defined in Section ^. H and computed in Q2UJI . For all x£l 2 , p(tq, n, x) 
has a limit p(0, ri,x) and, from Lemma El and Remark [2j we easily obtain the two- 
dimensional distributions of /(Co, ro, ri). 

Corollary 2. For all x, y £ M 2 and for all tq > 0, (/(Co,ro,ri,x),/(Co,ro,ri,y)) is a 
mixture of the two (two-dimensional) random variables (Co(x), Co(y)) and (Co(x), Ci(y)), 
where C\ is an independent copy of Co, with respective weights p(ro,ri,y — x) and 1 — 
p(ro,n,y - x). 

Note furthermore that, as in the isotropic case (|21[) . x t—> p(0, ri,x) is easily obtained 
by taking the limit as ro — > in the integrals of the RHS of (|2()[) : the one in the numerator 
is always finite by ()19|) and the one in the denominator is either positive or infinite by 
applying Lemma El More precisely, we have, for all x£K 2 , 

p(0,ri,x) = l(x = 0) if a > 3, (22) 
fn 1 7(x/tt) u 2 ~ a du , s 

^' ,)= Wo)-W.))^* " 6(U1, (23) 

By Corollary [21 (|2*2*)l implies that /(Co,0,ri) has independent (identically distributed if 
Co is stationary) samples if a > 3. For 1 < a < 3, (|23j) and Proposition El give that 
x 1 — > p(0, n, x) is a continuous M 2 — ► [0, 1] mapping. 

4.5 Preliminary properties of the limit process 

We have seen that, for a > 3, /(0, r±) is a white noise random field. On the contrary, for a € 
(2,3), Proposition ISIflij) shows that the bi-dimensional distributions given in Corollary El 
(taking the constant field for Co) exhibit interesting scaling properties. We have so far 
only been interested in finite-dimensional distribution of the colored dead leaves model. 
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Let us investigate how the simple scaling properties of the two-dimensional distributions 
influence the sample paths properties of the model. The first property we may check is 
the existence of a measurable version of I(Co,0, n). A random field {Z(uj,x) : x G R 2 } 
defined on (Q,G,P) is measurable if (w,x) h-> Z(w,x) is a (0 x R 2 , g®B(R 2 )) -> (R,B(R)) 
(jointly) measurable mapping (see e.g. |ST94| Section 9.4]). Recall also that the random 
field {Z(x), x G R 2 } is said to be stochastically continuous if, for all x G R 2 , Z(y) — Z(x) 
(Z(y) converges to Z(x) in probability) as y — > x. 

Proposition 12. Take a G (1,3) and r\ < oo. Assume that Co is stochastically con- 
tinuous. Then I (Co, 0,ri) is stochastically continuous. If moreover Co has a measurable 
version or, equivalently, if, for all G,H G B(R) and for all x G M 2 , y i— > Pr(Co(x) € 
G, Co(y) € i?) is a (R 2 , B(R 2 )) — > ([0, 1],B([0, 1])) measurable mapping, then there exists 
a measurable modification ofI(Co,0, r±), that is a measurable process /(Co, 0, n) defined on 
the same probability space such that, for all x G R 2 , Pr(I(Co, 0, ri, x) = /(Co, 0, n, x)) = 
1. 

proof For convenience we write I for /(Co,0,ri) in this proof. The two-dimensional 
distributions of / are given in Corollary |2j We have seen that x i— > p(0,ri,x) defined 
by Ij23(l is a continuous mapping. We obtain, for all x,y£l 2 and for all e > 0, 

Pr(|/(x)-/(y)| >e) <Pr(|C (x)-C (y)| >e) + (l-p(0,r 1 ,y-x)), 

which tends to zero as y — > x for Co stochastically continuous and since p(0, ri,0) = 1. 
Hence the first part of the proposition. For the second part, we apply [ST941 Theorem 
9.4.2]. We have to check two conditions on I, namely 

(i) There exists a countable set ScR 2 such that for all x£l 2 , there exists a sequence 
(x/j) such that /(x/%) — ► ^( x ) as k — > oo and x^sS for all k. 

(ii) For all G,H € B(R) and for all x G R 2 , y i-» Pr(/(x) G G, J(y) G if) is a 
(M 2 ,S(M 2 )) -» ([0, 1],S([0, 1])) measurable mapping. 

Condition (0) is a consequence of stochastic continuity. Moreover, condition (juj) must be 
satisfied by the random field Co if it has a measurable version (and since Co satisfies (0) 
as a consequence of stochastic continuity, condition (jnj on Co implies the existence of 
a measurable version for Co). Condition (JnJ) is easily checked on / by using the two- 
dimensional distributions given in Corollary [21 □ 

Corollary 3. The colored dead leave model 1(0, r±) is stochastically continuous and it 
admits a measurable modification if and only if a G (1,3). 
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proof For a > 3, 1(0, r\) is a white noise and thus does not have a measurable version 
(see |ST94| Example 9.4.3]) and neither is stochastically continuous. For a € (1,3) the 
constant random field trivially satisfies the assumptions of Proposition 1121 □ 

Let us conclude this section by noting that this result is coherent with a result of 
D. Mumford and B. Gidas which, in a paper investigating stochastic models for natural 
images, jMGOlj . show that a non-constant scale-invariant twodimensional random field 
cannot be a random function, and has to be a random distribution (generalized function). 
In the framework of the dead leaves model, this corresponds to the case a = 3, where there 
is no measurable version of the model. Indeed, the scale invariant case would correspond 
to a = 3 when ro — > and n — > oo simultaneously. In this case the limit of the two- 
dimensional distribution (/(x, ro,ri),/(y, ro,n)) depends on the ratio between logro and 
logri, but never on (x, y) (see Proposition IjnTjl) and we find a measurable limit only 
in the case where this limit is a constant field. In the next section, we investigate the 
regularity of our limit model in the framework of Besov spaces. 

Before proceeding, we illustrate the properties of the limit model with some simula- 
tions. In figure El we show two examples illustrating Proposition Images are simulated 
using the perfect simulation methods explained in Section |3J and gray levels are uniform 
between and 255. In the first example (left) we illustrate point (0); a = 2.5, and r\ — > oo. 
The image is of size 1000 x 1000, ro = 1, n = 100000; the process converges to a constant 
function. In the second example (right), we illustrate point (jnj); a = 3.5, and ro — > 0. The 
image is of size 10000 x 10000, ro = 1, r\ = 10000, the process converges to white noise. 
In figure ^ we illustrate the convergence of I(ro,ri) when ro — > and a = 2.9. The first 
image is of size 10000 x 10000, ro = 1, ri = 10000. The next three images are zooms of 
the same realization of the model (the zoom factor is two each time) . 

4.6 Some sample paths properties 

If ro > 0, the process /(Co, ro, r\) have paths for which occlusions influence the smoothness 
in a rather simple way: it introduces discontinuities along dM(ro,ri) and, between these 
discontinuities, the smoothness is driven by the properties of Co- For instance, in the 
case of a constant Co, J(ro,ri) simply is a piecewise constant process, the pieces being 
connected components of M?\dM(ro, ri). In the simple case where dY has finite length, 
dM(ro, ri) has locally finite length almost surely (as an easy consequence of Lemma@} and 
thus, under suitable regularity assumption for Co, is locally of bounded variation. This 
is contradicted by empirical experiments. In |GM01j . by investigating the distribution of 
sizes of "bilevel" sets in natural images (up to the smallest available scale), it is shown 
that this smoothness model is erroneous. In practice this means that a denoising approach 
relying on the a priori of a piecewise smooth image may interpret small objects as noise 
and, therefore, may result as a non-negligible loss of information. This is well known in 
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Figure 3: Illustration of the degenerate cases of Proposition [7| Left: case (0), the process 
converges to a constant (random) function. Right: case (jn)), the process converges to a 
white noise. 



the field of image processing. For instance, variational methods in the space of functions 
with bounded variation, such as the famous Rudin-Osher denoising scheme, |ROF92j . are 
known to erase textured area. The efficiency of these methods is dued to their ability 
to preserve the large scale geometric structure of images, but they often fail in handling 
the small scales structure of the image. A recent approach to overcome this difficulty has 
been proposed by Y. Meyer, see Mey02 , introducing a new functional space to account 
for textured regions in images, see also |OSV03| . [AABFC03] . 

The idea here is to propose an a priori model which includes the existence of small 
objects. How this model could be used efficiently in practice, say for denoising, may be 
a difficult question. In non-parametric regression (including the denoising problem), an 
important issue is to model the smoothness of the target function in terms of standard 
smoothness classes and, in many cases, this issue is closely related to an a priori statistical 
model (see e.g. |ASS98j ). The object of this section is to provide some insights in the 
smoothness properties of the model 1(0, r±) which we believe to be a promising a priori 
statistical model for natural images. We will consider Besov smoothness spaces, which have 
been standard spaces in nonparametric regression since it was promoted in the celebrated 
work of Donoho and Johnstone (see e.g. |DJ98j ). Of course not every value of a is relevant 
for the model 1(0, n). We will avoid a > 3 be cause in this case 1(0, ri) simply is a white 
noise (see Corollary I3J). In fact, the most suitable range for the modeling natural images by 
1(0, ri ) is 2 < a < 3. Indeed, in this case p(0,r±,-) exhibit interesting scaling properties 
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Figure 4: Illustration of the convergence of the model when a = 2.9 and tq — > 0. The first 
image (up left) is a realization when ro = 1, r\ = 10000, on a window of size 10000 x 10000. 
The three next images are zooms on this realization, the zoom factor being two each time. 
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at small scales (see Proposition [SJ which correspond to the empirical observations (see 
|Bud97| ) . Hence in the following we will only consider this case when ro = 0. When 
ro > 0, we will consider the usual case a G (1,3] and < ro < r\ < oo or a > 3 and 
< ro < r\ < oo. 

Since in all these cases, there exists a (M 2 ,£>(K 2 )) — > (R, Z3(IR)) measurable modifica- 
tion x i— ► I(lu,x) of I(ro,ri), it makes sense to investigate whether I belongs to some 
given functional spaces. We only take interest in local smoothness. Since the process is 
stationary, we may consider its restriction to the cube [0, l] 2 without loss of generality. 
We will take interest in Besov spaces and, as a byproduct, obtain that our limit process 
1(0, ri) is not of bounded variation, which is coherent with the results of |GM01j in view 
of natural images modeling. Let s G (0,1), p G [l,oo] and q G [1, oo]. The Besov space 
Bp' g ([0,l} 2 ) (see e.g. |DL9T] or |Mey90| ) is the Banach space endowed with the following 



norm 

i 

q d,U\l 



\f\B^--=\f\p + ( / Hf,u) p u- s y-\ 

\Ju>0 u J 

where is the usual L p ([0, l] 2 ) norm and to(f, u) p is the L p ([0, l] 2 ) modulus of smoothness 
of / at scale u, that is uj(f,u) p := sup| y | <n |A(/, y)| p , where A(/, y) is the difference 
operator applied to / with step y on [0, l] 2 , that is, the mapping x i— > (/(x + y) — 
/(x))l(xG[0,l] 2 ,x + y€[0,l] 2 ). 

Finally let us note that the results below can be easily generalized to a colored dead 
leaves process I(Cq, ro,r\) with Co satisfying the assumptions of Proposition 1121 but, in 
general, they would depend on Cq. Here we focus on the basic properties implied by the 
geometrical construction of the model. Therefore we only consider the case of 7(ro,ri), 
that is with Cq being the constant field, because then the variations of the colors are 
directly connected to the geometry of the model. For similar reasons, if BY has some 
irregularity, it may influence the smoothness of I{rQ,r\). This is clear from the example 
given in the end of Section 14.11 Hence the additional assumptions (A-EJ) or (A-I2J in the 
following results. 

Proposition 13. For all p G [1, oo] and for all s G (0,1), the following assertions hold 
true. 

(a) Let 1 < a < 3 and < ro < r\ < oo or if a > 3 and < ro < r\ < oo. Under (A\^), 



\T\ P 



E 

(b) Let 2 < a < 3 and = r < r\ < oo. Under (A^Bj), 

E 



< oo 4» s < 1/p. (24) 



\T\ P 



3 — a . . 

< oo s < . (25) 

p 
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proof In this proof we write A x B if there exists a constant c (possibly depending on 
the constants s, p, a,ro, r\ and 7) such that \B < A < cB. It is more convenient to use 
the modified modulus of smoothness 



«;(/,«)„:= (IV O / |A(/,y)| p dy. 

J\y\<u 

which satisfies w(-,-) p x oj{-,-) p (see |DL91j ). Because / is measurable, we may use the 
Fubini Theorem. We obtain 

E[\lf BS , P ] x E[\I\ p f + [ [ E[\A(I,yW p ](lVu- 2 )u-^dy- 

p Ju>oJ\y\<u u 

xl + / / (l-p(r ,r 1 ,y))(lVu- 2 )n^dy — 

Ju>oJ\y\<uAl u 

= 1 + 2 / / — — ; ; NN (IV u 2 )u ps dy — . (26) 

Assume we are in the case • By Lemma El the denominator and the numerator of 
the fraction in the RHS of Q26JI are continuous mapping of y G K 2 . Moreover, since 
< 7 < 7(0) the denominator is bounded away from zero and infinity independently of 
y. Under (A-J2J, the numerator x (1 A |y|). Case (jlj is then easily achieved. 

We now consider case (|b|). By Lemma El and since a < 3, the denominator writes 

1 (27(0) - %y/v)) v 2 ~ a dv x 7(0)77-7(3 - a) > 0. (27) 





Concerning the numerator, a change of variable gives, for all y / 0, 

r (7(0) - 7(y/«)) ^ a = |y| 3 ^ r /lyl ( 7 (0) - 7(y/(|y|t>))) ^ 
Jo Jo 

Consequently, since a < 3, under (A-EJ), 

/>oo 

sup / (7(0) — j(z/v)) v 2 ~ a dv < 00. 

|z|=l Jo 

By p9|). we have, for all z such that |z| = 1 and for all r £ (0, l/(2a2)), 



T(7(o) - 7(y/(|y W)) ^ 2 ~" d« = f 7(0) ^ 2 ~ Q rf« 

JO JO 



7(0) , 
3 — a 



,3-a 



By Lemma we know that 7(0) — 7(y/(|y|v)) is non-negative, so that the last three 
equations finally give 



(7(0) - i{y/v)) v 2 ~ a dv x |y| 3 "°(|y| V P a ~ 3 
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From (|26|). ()27|) and the last equations, we obtain 

E[\I\ P BS , P ] x 1 + / [ / |y| 3 " a dy^) (1 V n" 2 ^ -. 

p Ju>0 \J\y\<uM J u 

Hence the result. □ 

An almost sure smoothness result follows immediately. 

Corollary 4. For all p E [l,oo], for all q € [1, oo] and for all s £ (0,1), the following 
assertions hold true. 

(a) If 1 < a < 3 and < ro < ri < oo or if a > 3 and < ro < T\ < oo, under (A\Q). 

s <l/p^ I e B s p ' q a.s. (28) 

(b) 7/2 < a < 3 and = r < ri < oo, wider (04^, 

s < ^=-^ / G £^ a.s. (29) 

proof We simply use that if a non-negative random variable has finite expectation, it 
is necessary finite almost surely. Then well known inclusions of Besov spaces give the 
claimed results. □ 

It is worthwhile to observe that if / € Bp P in the parameter region s > l/p, then / 
is continuous. Hence the condition for Besov smoothness given in case (jc5$ is the almost 
(because the boundary s = l/p is not included) weakest condition for an almost surely 
discontinuous I. In other words the Besov smoothness condition (|24|) do not distinguish 
I with any other standard discontinuous random field from a Besov smoothness point 
of view. In particular it does not depend on a. In contrast the case ((b)) exhibits some 
interesting behavior. This behavior is of course closely related to the asymptotical ex- 
pansion found in Proposition |H1 in the isotropic case. Moreover in this interesting case, 
as we anticipated before, the range of Besov spaces in which lies 1(0, n) is strictly bigger 
than the space BV{[0,1} 2 ). Indeed, it is well known that BV{[0,1} 2 ) C B™' 1 C B s { 1 for 
any s > 1, and moreover that for any / G BV([0,1] 2 ), \\f\\ B v > C||/|| B oo,i > C"||/|| b m. 
Therefore, as a consequence of Proposition EI] ®> we have that, for any n > 0, a > 2, 
£||J(0,ri)|| BV = oo. 

Concluding remarks In this paper we introduced two dimensional random fields ob- 
tained as limits of sequences of dead leaves models. As advocated earlier, we believe them 
to be promizing a priori models for natural images. First it accounts for the geometry of 
natural images (in short the occlusion phenomenon). Second, we used a distribution of 
objects sizes based on empirical statistical properties of images, as promoted in previous 
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works . Third, getting rid of small sizes "cut off" allows for regularity results that are also 
compatible with empirical observations. 

Now, how this model can be used in practice remains open. On the one hand, the 
almost sure smoothness properties may be taken as a functional a priori for images. For 
instance, in a denoising framework, Corollary 0] provides a class of smoothness hypotheses 
only driven by parameter a. For more elaborate tasks such as shape extraction, it is clear 
that parameter a is not sufficient, and that geometrical properties of the model, depending 
on the distribution of objects shapes, should be investigated. To this end further work 
has to be done on the dead leaves model. On the other hand, the model could be used as 
a Bayes a priori, but this would need tractable knowledge of the dependance structure of 
the model. 
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